Accurate predictions for heavy quark jets 



Giulia Zanderighi(l) 
(1) CERN, Geneve 23, CH-1211, Switzerland 




Heavy-flavour jets enter many of today's collider studies, yet NLO predictions for these quan- 
tities are subject to large uncertainties, larger than the corresponding experimental errors. 
We propose a new, infrared safe definition of heavy-quark jets which allows one to reduce 
theoretical uncertainties by a factor of three. 



1 Introduction 

When looking at the current comparison between the inclusive 6-jet spectra measured by CDF 
and the corresponding next-to-leading order (NLO) predictions, Fig. [ipJ, one notices two striking 
features. Firstly, one sees a tension between data and theory: the ratio of data over NLO is 
around 1.2-1.5 over the whole range of accessible transverse momenta pt of the jets. Secondly, 
one notices that the uncertainties associated with the theoretical predictions are embarrassingly 
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Figure 1: Ratio of the measured inclusive &-jet spectrum to the NLO prediction. The measurement is performed 
for jets with transverse momentum 38 GeV < Prjet < 400 GeV and rapidity |yjet| < 0.7. 
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Figure 2: Top: X-factor for inclusive 6-jet spectrum as computed with MCFM, clustering particles into jets 
using the kt jet-algorithm with R=0.7, and selecting jets in the central rapidity region {\y\ < 0.7). Middle: 
scale dependence obtained by simultaneously varying the renormalisation and factorisation scales by a factor two 
around pt, the transverse momentum of the hardest jet in the event. Bottom: breakdown of the Herwig inclusive 
6-jet spectrum into the three major hard underlying channels contributions (for simplicity the small bb — > bb is 

not shown). 



large (~ 40—50%) for a NLO calculation and in particular they are larger than the corresponding 
experimental uncertainties. To understand why this happens it is useful to examine Fig. [21 

The top plots show that the large uncertainty is associated with very large K-factors. The 
middle plots confirm that the uncertainty is the same both with MCFlvPand MCNLd^. Finally, 
the bottom plots illustrate the origin of the poor convergence of the perturbative expansion: 
when breaking down the Herwig*^ 6-jet spectrum into the hard underlying channels it turns out 
that two NLO channels, flavour excitation, where a 6-quark is kicked out of the sea-quarks, and 
gluon splitting, where a gluon in the final state splits into a 66-pair, are larger than the leading 
order heavy quark production mechanism, flavour creation, when two incoming light partons 
produce a heavy quark pair. 

The reason why supposedly higher order contributions are actually larger than the leading 
order channel can be clarified by counting soft and collinear logarithms associated with the split- 
ting of gluons into 66-pairs. It turns out that flavour excitation contributes with (a^ Inp^/mfe)"' 
and gluon splitting contributes with (a^ Inp^/m;,)^""^ relative to the leading order, 0{a1) pro- 
cess. Since rrn, <C pt these contributions are enhanced. Moreover, the dominant contribution to 
the b-jet spectrum comes from jets originated from gluon splitting, which do not correspond to 
one's intuitive physical idea of a b-jet, one where a hard b is produced directly in the hard scat- 
tering. In the following we suggest to adopt a different jet-clustering algorithm to reconstruct 
6-jets. One that by making explicit use of the flavour information eliminates all higher-order log- 
arithmic enhancements associated to gluon splittings in the 6-jet spectra. This means that, after 
resumming initial state collinear logarithms into 6-pdfs, 6-jets can be computed using massless 
QCD calculations^ as long as one neglects power corrections m'^/pf (potentially log-enhanced). 



2 The heavy-quark jet algorithm 

We summarize here the inclusive heavy-flavour jet algorithm for hadron-hadron collisions'^. For 
any pair of final-state particles i, j define a class of distances d^J'"^ parametrized by < a < 2 
and a jet radius R 

^{F,a) _ ^ i ma.x{ku, hj)" mm{kti, ktj)'^~°' , softer of i,j flavoured, 
i?2 I min(A;|j, fc^ ) , softer of i,j flavourless. 



"We recall that according to the current experimental definition of a 6-jets, a b-jet is any jet containing at least 
one 6. 



where Rfj = Ayfj + A(t>fj, Ayij = Hi - Uj, A4>ij = 4>i - 4>j and ku, yi and (pi are respectively the 
transverse momentum, rapidity and azimuth of particle i, with respect to the beam. For each 
particle i define a distance with respect to the beam B at positive rapidity, 

^{F,a) ^ \ max{kti,hB{yi))°'m.m{ku,ktB{yi)y~°' , i is flavoured, 
*^ 1 mm{k^^, k^^{yi)) , i is flavourless, 

with 

ktBiy) = E i^ivi -y) + e{y - y^)ey^-'^) . (3) 

i 

Similarly define a distance to the beam B at negative rapidity by replacing ktB in eq. ([2]) with 

ktBiy) = E my - yi) + ®{y^ - y)^"'"') ■ (4) 

i 

Identify the smallest of the distance measures. If it is a d^^''^\ recombine i and j into a new 

particle, summing their flavours and 4-momenta; if it is a (or rf^^'"^) declare i to be a jet 

and remove it from the list of particles. Repeat the procedure until no particles are left. We 
define the 5-flavour or generally the heavy-flavour of a (pseudo)-particle or a jet as its net heavy 
flavour content, i. e. the total number of heavy quarks minus heavy anti-quarks. 

The IR-safety of this algorithm was proved in^. Apart from allowing one to take the limit 
rnq for the heavy quark mass (as long as collinear singularities associated with incoming 
heavy quarks are factorized into a heavy quark PDF), it ensures that one obtains the same 
results whether one considers heavy-quark flavour at parton level, or heavy-meson flavour at 
hadron level, modulo corrections suppressed by powers of ^qcd/pi- 



3 Results 

Our results are summarized in fig. [sP where we show the inclusive 6-jet pf-spectrum as obtained 
with the flavour algorithm specified above with a = 1, and R = 0.7, the latter having been 
shown to limit corrections associated with the non-perturbative underlying event The left 
(right) column of the figure shows results for the Tevatron Run II (LHC). We have selected 
only those jets with rapidity |y| < 0.7. We also show the full inclusive jet spectrum (all jets) as 
obtained with a standard inclusive /ct-algorithm^S] with R = 0.7. 

We notice the considerable reduction of i^-factors, which are around 1.3 and the moderate 
uncertainties associated with scale variation, signaling that the perturbative expansion is now 
well under control. Our predictions constitute therefore the first accurate predictions for inclusive 
heavy quark jets. 

We remark that very similar results are obtained when considering charmed jet spectra. An 
interesting issue there is that predictions are very sensitive to possible intrinsic charm compo- 
nents of the proton^^. This means that this type of observable has a potential to set constraints 
on such intrinsic components. 

A last remark concerns the feasibility of the experimental measurement of heavy flavour 
jets defined with our flavour algorithm. Our jet-clustering algorithm requires that one identify 
heavy-fiavoured particles and that one uses a different distance measure when clustering heavy 
or light objects according to eq. ([1]). It is particularly important to identify cases when both 
heavy flavoured particles are in the same jet, so as to label this jet a gluon jet and eliminate it 
from the 6-jet spectrum. Experimentally techniques for double 6-tagging in the same jet already 
exist and steady progress is to be expected in the near future^^^T However one has always a 
limited efficiency for single h tagging, and even more for double 6-tagging in the same jet. On 
the other hand preliminary studies indicate that one does not necessarily need high efficiencies. 
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Figure 3: Inclusive jet spectrum at the Tevatron (right) and at the LHC (left). The top two panels sho-w results 
for both 6-jets and all-jets, -while the lo-wer three panels apply only to 6-jets. See text for further details. 



but -what is more crucial is that one understand those efficiencies wehf^. We look for-ward to 
further investigation in this direction. 
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